Stochastic Acceleration of ^He and ^He by Parallel Propagating 

Plasma Waves 

Siming Liu^, Vahe Petrosian^'^, and Glenn M. Mason'^ 
ABSTRACT 

O 

o 

CN ' Stochastic acceleration of ^He and ^He from a thermal background by parallel 

■ propagating turbulent plasma waves with a single power-law spectrum of the 

pLn ' wavenumber is studied. In the model, both ions interact with several resonant 

CX5 ' modes. When one of these modes dominates, the acceleration rate is reduced 

^ ■ considerably. At low energies, this happens for ^He, but not for ^He where 
contributions from the two stronger modes are comparable so that acceleration of 
^He is very efficient. As a result, the acceleration of ^He is suppressed by a barrier 

O ■ below ~ 100 keV nucleon"^ and there is a prominent quasi-thermal component 

ff^ ■ in the "^He spectra, while almost all the injected ^He ions are accelerated to 

^ ' high energies. This accounts for the large enrichment of ^He at high energies 

O ■ observed in impulsive solar energetic particle events. With reasonable plasma 

^ ■ parameters this also provides a good fit to the spectra of both ions. Beyond 



> 



1 MeV nucleon ^, the spectrum of ^He is softer than that of ^He, which is 



2 ■ consistent with the observed decrease of the ^He to ^He ratio with energy. This 

^ ■ study also indicates that the acceleration. Coulomb losses and diffusive escape 

of the particles from the acceleration site all play important roles in shaping the 



ion spectra. This can explain the varied spectral shapes observed recently by the 
Advanced Composition Explorer. 
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1. Introduction 

The abundance and spectra of most ions observed in impulsive solar energetic particle 
events (SEPs) qualitatively agree with predictions of stochastic acceleration (SA) by plasma 
waves or turbulence (PWT) (e.g. Mobius et al. 1980; 1982). However, the observed enrich- 
ment of ^He, sometimes by more than four orders of magnitude over its coronal abundance 
(Hsieh & Simpson 1970; Serlemitsos & Balasubrahmanyan 1975; Mason et al. 2002), re- 
mains a theoretical challenge. The aim of this paper is to provide an explanation for this 
observation in the frame of SA models. 

Solar ^He rich events are enriched also in other heavy elements (Hurford et al. 1975; 
Mason et al. 1986) and are closely related to impulsive X-ray flares, scatter-free electron 
events, and type III radio bursts (Reames, von Rosenvinge & Lin 1985; Reames, Meyer & 
von Rosenvinge 1994; Mazur, Mason & Klecker 1995; Mason et al. 2002), indicating that 
a common particle acceleration process may be responsible for all these features. Among 
the acceleration mechanisms applied to solar flares, the theory of SA, a second order Fermi 
acceleration process by PWT, has achieved several signiflcant successes. The theory was 
flrst introduced to account for the acceleration of protons and other ions responsible for 
the observed nuclear gamma ray line emissions (Ramaty 1979; Hua, Ramaty & Lingenfelter 
1989). Later it was shown that it is also apphcable to the acceleration of electrons, which 
produce the impulsive hard X-ray emission (Miller & Ramaty 1987; Hamilton & Petrosian 
1992; Park, Petrosian & Schwartz 1997; Petrosian & Donaghy 1999). A recent study of SA 
by cascading Alfvcn turbulence indicated that it may also explain heavy ion acceleration as 
well (Miller 2003). 

It is generally accepted that the enrichment of ^He is related to its peculiar charge to 
mass ratio so that certain plasma waves resonantly accelerate or heat it to high energies, 
producing the observed abundance ratio at a few MeV nucleon"^ (Mobius et al. 1980; 
Zhang 1999). Previous work either discusses the feasibility of selective accelerating ^He 
from very low energies (Temerin & Roth 1992; Miller & Vihas 1993), or considers a two- 
stage mechanism; a selective preheating process followed by an acceleration applicable to 
all ions (Fisk 1978). The former cannot be compared with observation directly because the 
acceleration of other ions has not been addressed accordingly. For the latter, the observed 
abundance pattern at high energies is determined by the preheating processes (Paesold, 
Kartavykh & Benz 2003) and the predicted spectrum depends on several parameters with 
the corresponding physical processes unspecifled (Sakurai 1974; Zhang 1995). 

Given its simple dispersion relation, Alfven turbulence is used in most of the SA models 
(Ramaty 1979; Mobius et al. 1982), which cannot explain the varied ion spectral forms and 
underestimates the production rate of higher energy particles in disagreement with the recent 
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observations by the Advanced Composition Explorer {ACE; Mason et al. 2002). However, low 
energy particles mostly interact with cyclotron and/or whistler waves. To understand the 
acceleration of particles from a thermal background, it is necessary to use the exact dispersion 
relation (Stcinacker et al. 1997). Wc have carried out such an investigation recently, which 
addresses the acceleration of electrons vs protons by waves propagating parallel to magnetic 
field lines (Petrosian & Liu 04, PL04 for short). Besides its achievements in explaining 
many features of solar flares, we found that the acceleration of low energy protons can be 
suppressed significantly relative to the electron acceleration due to the dominance of the 
wave-proton interaction by the resonant Alfven mode. 

In this letter, we show that a similar mechanism makes the acceleration of low energy 
^He much more favorable relative to that of ^He. Because their loss rates are comparable, 
such a mechanism basically depletes ^He from the thermal background to ~ MeV energies 
while most of the ^He ions are "trapped" at low energies, resulting a prominent quasi-thermal 
^He component. The model not only explains the enrichment of ^He at a few MeV nucleon"^, 
but also accounts for the flattening of the ^He spectrum at lower energies (Mason, Dwyer 
& Mazur 2000; Mason et al. 2002). It also predicts that beyond a few MeV nucleon"^ the 
spectrum of ^He is always harder than the ^He spectrum. In the next section, we present 
our SA model. Its application to SEPs is discussed in § 3. In § 4, we summarize our results 
and discuss their implications. 

2. Stochastic Ion Acceleration 

The theory of SA by parallel propagating waves is explored in detail in PL04. Here, we 
highlight a few key points, which are crucial in understanding the ion acceleration process. 

First, a magnetized plasma can be characterized by the ratio of the electron plasma 
frequency a;pe to gyrofrequency Qe- 

a = ojpjn^ = 3.2(ne/10^°cm-3)i/2(5o/100 G)"^ , (1) 

where Ue is the electron number density and Bq is the large scale magnetic field. To describe 
the wave modes in the plasma, one also needs to know the abundance of ions. In our case, 
the inclusion of a-particle, whose number density is about 10% of the proton density, is 
essential while other heavier elements have little effect on the dispersion relation ui = uj{k) 
(Steinacker et al. 1997), where lo is the wave frequency and k is the wavenumber. 

Low energy ions mostly interact with left-handed polarized waves. Figure 1 shows the 
dispersion relation of these wave modes as indicated by the dotted curves in a plasma with 
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a — 0.45 and Yne = 0.08, where Yne is the abundance of ct-particle (here following the nota- 
tion in PL04, negative frequencies indicate that the waves are left-handed polarized). There 
are two distinct branches, which we call ^He-cyclotron branch (HeC: upper one) and proton- 
cyclotron branch (PC: lower one) because they asymptotically approach the corresponding 
ion cyclotron waves at large wavcnumbers. At small A;, HcC branch gives the Alfven waves, 
while the PC branch approaches a wave frequency cupc — (1/2 + YHe)^^p, where f2p is the 
proton gyrofrequency. 

The acceleration process is determined by the resonant wave-particle interaction. The 
waves and particles couple strongly when the resonant condition uo = kjScfi — Qi'~f~'^ , is 
satisfied, where i3c, ji, and 7 are respectively the velocity, the pitch angle cosine, the 
nonrelativistic gyrofrequency, and the Lorentz factor of the ions. Figure 1 also shows this 
resonant condition (straight lines) for ^He and ^He with ji — 1 and two values of energy: 
E — 0.5 and 47 keV per nucleon. The intersection points of these fines with the dotted curves 
for the dispersion relation satisfy the resonant condition and indicate strong wave-particle 
couplings. Tfie wave-particle interaction rates are calculated by adding contributions from 
these points (PL04). 

At low energies (e.g. E — 0.5 keV nucleon"^ in the figure) ^He interacts only with two 
waves, one from the HeC and one from PC (not shown in the figure because of its large 
wavenumber). On the other hand, '^Hc can interact with four waves, three of them from 
the PC branch (two of them are shown in the figure) and one from the HcC branch (also 
not shown in the figure)^. As shown in Figure 1 there arc two resonant modes for '^Hc with 
nearly equal and small wavenumbers while "^He has one resonant mode at low k (see circles 
in Fig. 1); the next important resonant mode of ^He is at much larger k value beyond the 
range of the figure. Consequently, for a power law spectrum of turbulence, the interaction 
of low energy ^He will be dominated by a single mode from the HeC branch, giving rise 
to a very low acceleration rate relative to ^He (see discussion below) ^. This is not true for 
particles with high energies [E — 47 keV nucleon"^ in the figure). For these particles, the 
two stronger resonant modes have comparable contributions for both ions. 

Following previous studies (Dung & Petrosian 1994; PL04), we will assume that the 
turbulence has a power law spectrum £{k) oc A;~^ with a low wavenumber cutoff /Cmin, which 



^Both ions also interact with a right-handed polarized wave in the electron cyclotron branch. Because its 
contribution to the interaction rates of low energy ions is small, we do not discuss this branch here. However, 
the contribution of this wave branch is included in the numerical calculation below. 

^It should be noted that this in not true for ions with nearly 90° pitch angle (/U <C 1). However, this 
contributes little to the acceleration as a whole because there are few such particles 
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corresponds to the turbulence generation scale and determines the maximum energy of the 
accelerated particles. For a turbulence spectral index g < 1, a high wavenumber cutoff is also 
required to ensure the convergence of the total turbulence energy density. In this case, we 
choose k^nax large enough so that the acceleration of particles from the thermal background 
is not affected by this cutoff. Then the characteristic interaction time scale Tp is given by 
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where £tot is the total turbulence energy density. 
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Ions in a turbulent plasma can gain energy by interacting with the waves with diffu- 
sion coefficients Dpp, D^p, and -D^^, where p is the momentum of the particles, and lose 
energy via Coulomb collisions with electrons and protons in the thermal background. As 
a result the particles diffuse in both real and momentum space. The spatial diffusion can 
be approximated with an escape term with a time scale Tcsc- Because the Coulomb scat- 
tering rate at low energies and the pitch angle scattering rate D^^ at high energies are 
larger than the acceleration rate (except for highly magnetized plasmas with a <C 0.02), 
the particle distribution is nearly isotropic over all energies. The energy diffusion rate 
l{l^:P) — '^Z^it^^P) — {Dpp — Dnp/Di^ij)/p'^ (see Schlickeiser 1989; Dung & Petrosian 1994). 
As stressed in PL04, when the interaction is dominated by one of the resonant modes, this 
rate becomes negligibly small and the particle acceleration is suppressed. In PL04, we showed 
that such a suppression occurs for low energy protons but not for electrons. As mentioned 
above at low energies the '^He interaction is dominated by one mode resulting in a low accel- 
eration rate relative to '^Hc. For isotropic pitch angle distribution, the particle distribution 
integrated over the acceleration region, N[e), as a function of the total (not per nucleon) 
kinetic energy £, satisfies the well known diffusion-convection equation: 

dN d"^ d N 

where D^^ = < ^{Hjp) >~ eA{e) is the pitch angle averaged diffusion rate, Q is a source 
term, and the loss rate e^L is given in PL04. 

Figure 2 shows the relevant time scales for a model with L = 10^^ cm, Bq = 500 
Gauss, -n-e = 5 X 10^ cm~^ and temperature T = 0.2 keV. The turbulence spectral index 
q = 1, = 3 s^^ and a = 0.45. The acceleration and loss times are defined as: Ta = 
e/A{s) ,Tioss = Depending on L, Tp and s, either the transit time Ttr = {L/\/2/3c) 

or the diffusion time rjis = 2t^^/tsc can dominate the escaping process, where Tgc is the 
scattering time including both wave-particle scatterings and particle-particle collisions. We 
define Tesc = {L/\/2(3c){l -\- \/2L/ (3cTsc) , which incorporates both processes of escape. 
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There are several features in these times, which require particular attention: 

First, all the time scales are energy dependent and the acceleration times have compli- 
cated forms, in contrast to the often used simplified acceleration by pure Alfven turbulence 
and/or with energy independent escape time. This affects the spectral shapes significantly. 

Second, the typical time scale varies from less than one second to hundreds of seconds 
with different processes dominating at different energies. 

Third, the acceleration times for "^He and ^He are quite different due to interaction 
with different waves, while their escape and loss times are similar because these times are 
dominated by Coulomb collisions and wave-particle scatterings (oc< >). 

Fourth, the acceleration time of '^He is shorter than the Coulomb loss time below a few 
hundreds of keV micleon"^. Consequently, for a background plasma of a few million Kelvin 
almost all the injected ^He can be accelerated to high energies (a few MeV nucleon"^). Above 
this energy the escape time becomes shorter than the acceleration time and the acceleration 
is quenched. 

Fifth, ^He acceleration time increases sharply with the decrease of energy near 100 kcV 
nucleon"^. This is because at lower energies the wave-particle interaction is dominated by 
one of the resonant modes. ^He acceleration time starts to decrease with energy near 50 keV 
nucleon"^ where the two resonant waves have comparable wavenumber (Figure 1). Because 
the acceleration time is longer than the loss time at low energies, for an injection plasma 
with a temperature less than ~ 10^ K, only a small fraction of ^He in the Maxwellian tail can 
be accelerated to very high energies. The rest is heated up into a quasi-thermal distribution. 

Sixth, the acceleration time of ^He is shorter than that of ^He above ~ 100 keV 
nucleon"^, so that the ^He acceleration becomes more efficient at these energies. 

3. Application to SEPs 

Because the duration of SEPs are much longer than the characteristic time scales dis- 
cussed in the previous section, we will assume that the system reaches a steady state. One 
can then ignore the time derivative term in equation (3) and solve this equation for the 
escaping flux / = N/T^ci which can be directly compared with the observed fluences. 

Figure 3 shows our model fit to the spectra of an SEP event on January 6, 2000, the event 
with the largest ^He enrichment (Mason et al. 2000). The data are obtained by carefully 
subtracting the background and avoiding contaminations from nearby events. For the model 
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spectra, we assume that the injected plasma has a coronal abundance, i.e. the abundance of 

■^He is two thousand times higher than that of ^He. The model parameters are the same as 
those in Figure 2. The temperature of 0.2 keV is a typical value for impulsive SEP events 
(Reames et al. 1994). Because the acceleration region is connected to open magnetic field 
lines probably above the flaring coronal loops, a size of L = 10^°cm seems quite reasonable. 
The magnetic field and gas density are also characteristic of an upper coronal reconnection 
region. The turbulence parameters q and Tp are free model parameters. 

The model not only accounts for the enrichment of '^He, but also gives a reasonable fit 
to the spectra. Because at low energies the acceleration time of ^He is shorter than any other 
times, almost all of the injected ^He ions are accelerated to a few hundreds of keV nucleon^^, 
where Ta becomes longer than the loss and escape time, resulting in a sharp cutoff. The ^He 
acceleration time, because of the barrier mentioned above, is much longer than its loss time 
below a few keV nucleon"^. This suppresses the ^He acceleration so that most of the injected 
^He ions remain near the injection energy and form a quasi-thermal component. However, 
even when A{e) < and the direct acceleration is difficult, some ^He ions can still diffuse 
to high energies due to the first term on the right hand side of equation (3). Consequently, 
we expect a nonthermal tail in addition to the quasi-thermal component. This tail is cut 
off at ~ 10 MeV nucleon"^ (Fig- 3), when the escape time becomes much shorter than the 
acceleration time (Fig. 2). The acceleration barrier extends to ~ 100 keV nucleon"^, most 
of the ^He ions diffuse to higher energies are stopped below this energy. About 10% of the 
particles can be accelerated to even higher energies and produce the observed high energy 
spectrum with some spectral features in this energy range. However, a more realistic model, 
including acceleration by other wave modes, is expected to produce smoother spectra. 

The above behaviors depend on six model parameters, which are shown in Figure 2 
(note that a is given by rig and Bq via eq.[l]). Figure 4 shows how the spectra change with 
two of these parameters: (or the energy density of the turbulence) and the injection 
temperature T. Because of the acceleration barrier at low energies, the acceleration of "^Hc 
is very sensitive to T. In the figure, we adjust Tp in order to make the fiucnce of '^Hc and 
^He comparable at a few hundred keV per nucleon. As evident, ^Hc spectrum peaks at a 
higher energy with the increase of the intensity of the turbulence and the ^He enrichment 
also changes significantly. This can explain the variation of the ^He spectral peak and its 
enrichment in class 2 SEPs (Mason et al. 2000; Mason, Mazur & Dwyer 2002). 
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4. Summciry and Discussion 

In the paper, we have shown that the observed "^He and "^He spectra in some impulsive 
SEPs can be produced via SA by PWT propagating parallel to magnetic field lines. In the 
model, the particles are injected from a thermal background. Because the interaction of low 
energy ^He is dominated by one of the resonant wave modes, an acceleration barrier forms, 
giving rise to a prominent quasi-thermal component in the steady state particle spectrum. 
This is not true for ^He ions, however, which can be efficiently accelerated by two waves in 
the PC branch, from very low energies to ~ MeV nucleon"^, where they escape quickly from 
the acceleration site, producing a sharp spectral cutoff. This explains the ^He enrichment 
and its varied spectrum seen in class 2 SEPs, which have relatively short rise-time and are 
more likely associated with the impulsive phase particle acceleration directly. Class 1 events 
have less ^He enrichment and similar spectra for most ions, indicating a different acceleration 
process (Mason et al. 2000). The feasibihty of producing such events by shocks or turbulence 
with a different spectrum than that of the present model will be investigated in the future. 

Because the "^He acceleration time is shorter than that of ^He at higher energies and 
their loss and escape times are comparable, the model predicts a decrease of the ^He to ''He 
ratio with energy beyond ~ 1 MeV nucleon"-*^, which is consistent with observations (Mason 
et al. 2002; Reames et al. 1997; Mobius et al. 1982; Mobius et al. 1980; Hsieh & Simpson 
1970). Spectral information at higher energies is critical in testing this model prediction. We 
emphasize here that a special injection mechanism is not required to explain the observations. 
However, a second-phase shock acceleration (Van HoUebeke, McDonald & Meyer 1990) can 
modify the ion spectra and may explain events where the high energy spectral indices of ^He 
and ^He are similar (Serlemitsos & Balasubrahmanyan 1975; Reames et al. 1997). 

The model can also be applied to acceleration of heavy-ions. In this case, the acceleration 
is dominated by resonant interaction with waves in the HeC branch and ions with low 
charge to mass ratio are accelerated more efficiently, which may explain the increase of ion 
enrichment with the increase of their charge to mass ratio (Mason et al. 1986). Because 
^He acceleration is dominated by interaction with the PC branch, which may be decoupled 
from the HeC branch, the model may also explain the lack of quantitative correlation of the 
enrichment of ^He and heavy- ions (Mason et al. 1986). We will investigate these aspects in 
future work. 

The primary uncertainties of the model are associated with the turbulence generation 
mechanism and the subsequent cascade and dissipation processes. We start the investigation 
with a power law spectrum of the wavenumber for the turbulence, which can be characterized 
by the turbulence energy density and the spectral index. The acceleration of particles by 
obliquely propagating waves and turbulence with complex spectral shapes may give quanti- 



-9- 



tatively different results. However, tfie idea of an efficient depletion of low energy ^He from 
the source region to high energies to produce the observed enhancement seems well founded. 

This research was partially supported by NSF grant ATM-0312344, NASA grants NAG5- 
12111, NAGS 11918-1 at Stanford and NASA grant PC 251429 at University of Maryland. 
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Fig. 1. — Dotted curves: dispersion relation of PC and HeC left hand polarized waves 
parallel to the large scale magnetic field. Solid lines: resonant condition for 47 keV/nucleon 
^He {thick) and ^He {thin). Dashed lines: resonant condition for 0.5 keV/nucleon. Circles 
designate the resonant points. 
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Fig. 2. — Time scales of a SA acceleration model. The model parameters are indicated in 
the figure. Solid: ^He, dashed: ^He. The thicker curves indicate the acceleration times. See 
text for details. 
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Fig. 3. — Model fit to the spectra of tlie SEP on Jan. 6, 2000. Tlie tliick solid curve is for '^He 
and the dashed curve is for ^He. The injected particles have a temperature of T = 0.2 keV 
and the corresponding distribution of ^He is indicated by the thin solid curve with arbitrary 
normalization. All other model parameters are shown in Figure 2. 
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Fig. 4. — Dependence of the ^He (thick curves) and ^He (thin curves) spectrum on Tp and 
T. The dotted, sohd, and dashed curves are for = 1, 3, and 8 s~^ and T = 1.7, 0.2, and 
0.09 keV, respectively. The dotted (dashed) curves are shifted up (down) by decade. All 
other model parameters remains the same as those in Figure 2. 



